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\ Abstract 

Q ■ We propose to describe the process pp — )• D°D° in a perturbative QCD motivated 

framework where a double-handbag hard process udud — > cc factorizes from transition 
^ \ distribution amplitudes, which are quasi forward hadronic matrix elements of ^q^^c op- 

erators, where q denotes light quarks and c denotes the heavy quark. We advocate that 
the charm-quark mass acts as the large scale allowing this factorization. We calculate this 
Q-[ process in the simplified framework of the scalar diquark model and present the expected 

qj cross sections for the PANDA experiment at GSI-FAIR. 



(N 



X 

5-H 



PACS: 13.75.-n, 25.43.+t, 13.85.Fb 



m : 1 Introduction 

Ov 

qq ■ The collinear factorization framework allows to calculate a number of hard exclusive ampli- 
tudes in terms of perturbatively calculable coefficient functions and non-perturbative hadronic 
matrix elements of light-cone operators. The prime example is the calculation of the deeply 
virtual Compton- scattering (DVCS) amplitude in the handbag approximation with generalized 
parton distributions (GPDs), non-forward matrix elements of a quark-antiquark non local op- 
erator v P(z) v P(0) between an incoming and an outgoing baryon state. Strictly speaking, this 
description is only valid in a restricted kinematical region, called the generalized Bjorken scal- 
ing region, for a few specific reactions, and in the leading-twist approximation. It is, however, 
suggestive to extend this framework to the description of other reactions where the presence of 
a hard scale seems to justify the factorization of a short-distance dominated partonic subprocess 
from long-distance hadronic matrix elements. Such an extension has, in particular, been pro- 
posed in [1] for the reaction pp — >■ A C A C with nucleon to charmed baryon GPDs. The extension 
of the collinear factorization framework to the backward region of DVCS and deep exclusive 
meson production [2, 3] leads to the definition of transition distribution amplitudes (TDAs) as 
non-forward matrix elements of a three-quark non-local operator v P(z) v P(v) v P(0) between an 
incoming and an outgoing state carrying a different baryon number. Here too, this description 
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is likely to be valid in a restricted kinematical region, for a few specific reactions, and in the 
leading twist approximation. We propose here to extend the approach of [1] to the reaction 
pp -> D°D° which will be measured with the PANDA [4] detector at GSI-FAIR. For this pro- 
cess the baryon number exchanged in the t-channel implies that hadronic matrix elements with 
v P(z) x P(y) l P(0) operators enter the game. Let us stress that we have no proof of the validity of 
this approximation but take it as an assumption to be confronted with experimental data. For 
this approach to be testable, one needs to model the occurring nucleon to charmed meson TDAs. 
In contrast to the N — > % TDAs, which have been much discussed [5], we do not have any soft 
meson limit to normalize these TDAs. We will rather use an overlap representation in the spirit 
of [6]. 



2 Hadron Kinematics 




Figure 1: Kinematics of pp — > D°D° in the symmetric CMS. 

The kinematical situation for pp — >■ D°D° scattering is sketched in Fig. 1. The momenta 
and helicities of the incoming proton and antiproton are denoted by p, jl and q, v, the momenta 
of the outgoing D° and D° by p' and q', respectively. The mass of the proton is denoted by 
m, that of the D° by M. We choose a symmetric center-of-momentum system (CMS) in which 
the longitudinal direction is defined by the average momentum of the incoming proton and the 
outgoing D°, respectively. The transverse momentum transfer is symmetrically shared between 
the incoming and outgoing hadrons. 
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In light-cone (LC) coordinates the hadronic momenta are parameterized as follows, 
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where we have introduced sums and differences of the hadron momenta: 

P-=\{p + p'), q-=\{q + q) and A := p' - p = q- q' . (2) 

The minus momentum components can be obtained by using the on-mass shell conditions p 2 = 
q 2 = m 2 and p' 2 = q' 2 = M 2 . The skewness parameter ^ gives the relative momentum transfer 
in the plus direction, i.e. 

p + + p + 2p + 

The Mandelstam variable s is given by 

s = (p + q) 2 = (p' + qf. (4) 

In order to produce a D°D° pair, s must be larger than AM 2 . The remaining Mandelstam vari- 
ables, t and u, read: 

t= A 2 = ( P '-p) 2 = (q-qf (5) 

and 

u = (q'-p) = (p - q) , (6) 

so that s + 1 + u = 2M 2 + 2m 2 . For later convenience we also introduce the abbreviations 

A m := ^\-Am 2 /s and A M := \J\-AM 2 js. (7) 
For further relations between the kinematical quantities see App. A. 



3 Double Handbag Mechanism (DHM) 

The double handbag mechanism (DHM) which we use to describe pp — > D°D° is shown in 
Fig. 2. It is understood that the proton emits an S[w<i]-diquark with momentum k\ and the anti- 
proton a S[w<i]-diquark with momentum ki. They undergo a scattering with each other, i.e. they 
annihilate in our case into a gluon which subsequently decays into the heavy cc pair. Those 
produced heavy partons, characterized by k\, X[ and k' 2 , Aj, are re-absorbed by the remnants of 
the proton and the antiproton to form the D° and the D°, respectively. One could, of course, also 
think of vector-diquark configurations in the proton and V [ud] V [ud] annihilation to produce the 
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Figure 2: The handbag contribution to the process pp — > D°D°. The momenta and helicities of 
the baryons and quarks are specified. 

cc pair. But in common diquark models of the proton it is usually assumed that the probability to 
find a V [ud] diquark is smaller than the one for the S [ud] diquark. Further suppression of V [ud] 
diquarks as compared to S [ud] diquarks occurs in hard processes via diquark form factors at 
the diquark-gluon vertices [7]. We thus expect that our final estimate of the D°D° cross section 
will not be drastically altered by the inclusion of vector-diquark contributions and we stick to 
the simpler scalar diquark model. 

The whole hadronic four-momentum transfer A is also exchanged between the active partons 
in the partonic subprocess 

S[ud](k l )S\ud](k 2 ) ->• c(k[, X[) c(k' 2 , l 2 ) . (8) 

In (8) we neglect the mass of the S[ud]— (anti)diquark, but take into account the heavy (anti-) 
charm-quark mass m c . In order to produce the heavy cc pair, the Mandelstam variable s of the 
partonic subprocess has to be 

s > Am 2 c , (9) 

where Am 2 c ~ 6.5 GeV 2 . We have taken the (central) value for the charm-quark mass from the 
Particle Data Group [8], which gives m c = 1.275 ± 0.025 GeV. Thus the heavy-quark mass 
m c is a natural intrinsic hard scale which demands that the intermediate gluon has to be highly 
virtual. This allows us to treat the partonic subprocess perturbatively, even at small — t, by 
evaluating the corresponding Feynman diagram. All the other non-active partons inside the 
parent hadrons are unaffected by the hard scattering and thus act as spectators. 
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For the double handbag mechanism the hadronic pp — > D°D° amplitude can be written as 

d A z\ * 171 f , 4r n/T _. f d 4 z 2 



(') a' 



(2k) 4 ~ J " ^"^> J (2k) a 
x(&:p'\£?X , a ,(-z 1 /2)4>f!" rfl (+z 1 /2)|p:p, At ) J y <»Jh,h) (10) 



x (D° : ?'|^^ lt (+z 2 /2)^ a ,(-z 2 /2)|p : q,v) 



where the assignment of momenta, helicities, etc. can be seen in Fig. 2. af' and a\ denote 
color and spinor indices, respectively. In analogy to the hadronic level we have introduced the 
average partonic momenta k\ := (ki + ty /2, i = 1,2, of the active partons. We note once more 
that the full hadronic momentum transfer is also transferred between the active partons, i.e. 
k[—k[ = p — p' = k' 2 — &2 = q' — q- 

The hard scattering kernel, denoted by H w ,(ki^2), describes the hard S[ud]S[ud] — > cc sub- 

process. The soft part of the p — > D° transition is encoded in the Fourier transform of a hadronic 
matrix element which is a time-ordered, bilocal product of a quark and a diquark field operator: 

/^^ lZl (^:p1^^ W{ (-^i/2)^ 1 (+zi/2)]p:p, J u>. dD 

In (11) 4> 5 [" d '(+zi/2) takes out an 5[«J]-diquark from the proton state \p : p,fl) at the space- 
time point zi/2. TheSlud] diquark then takes part in the hard partonic subproces. The x ¥ c (—zi/2) 
reinserts the c-quark at — zi/2 into the remnant of the proton which gives the desired final 
hadronic D° state \D° : p'). At this stage the appropriate time-ordering of the quark field opera- 
tors (denoted by the symbol ST) has to be taken into account. The remnant of the proton, which 
does not participate in the hard partonic subprocess, constitutes the spectator system. For the 
p — >■ D° transition we have the Fourier transform 

/ ^ /k2Z2{D ° : J\^*% d] \+Z2/2)% 2a ,(-Z2/2)\p :q,v), (12) 

which can be interpreted in an analogous way as (11). The pp — > D°D° amplitude (10) is thus a 
convolution of a hard scattering kernel with hadronic matrix elements Fourier transformed w.r.t. 
the average momenta k\ and k.2 of the active partons. 

For the active partons we can now introduce the momentum fractions 

x, := and x\ := -\ . (13) 

p+ p'+ 

For later convenience we also introduce the average fraction 

ki ~~\~~ k-i k~i 
p+j r pf+ p + 
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which is related to x\ and x[ by 



xi = 



1+e 



and 



Xi = 



(15) 



respectively. 

As for the processes in Refs. [1] and [9], due to the large intrinsic scale given by the heavy 
quark mass m c , the transverse and minus (plus) components of the active (anti-)parton momenta 
in the hard scattering kernel H are small as compared to their plus (minus) components. Thus 
the parton momenta can be replaced by vectors lying in the scattering plane formed by the parent 
hadron momenta. For this assertion one only has to make the physically plausible assumptions 
that the momenta are almost on mass-shell and that their intrinsic transverse components (di- 
vided by the respective momentum fractions (15)) are smaller than a typical hadronic scale of 
the order of 1 GeV. We thus make the following replacements: 



k\ ->■ 
k' 2 



r 2 A 2 



-xv 



8A, ' '" 2 



with A:^ =^i/? + , 
with 



with A: 2 =x 2 q 



m 2 c +4±\/4 _ , A 



2^ 



with fc 2 = jc 2 ^' 



(16) 



As a consequence of these replacements it is then possible to explicitely perform the integrations 
over k x , k 2 , kjj and k^ 2 - Furthermore the relative distance between the (anti-)S[w<i]-diquark 
and the (anti-)c-quark field operators in the hadronic matrix elements is forced to be light-like, 
i.e. they have to lie on the light cone and thus the time ordering of the field operators can be 
dropped. After these simplifications one arrives at the following expression for the pp —?■ D°D° 
amplitude: 



£ J dkte(k+) j ^ e fk ^ J dk 2 e(k 2 ) J 



/XV 



' dz 2 fir z+ 



2n J L v L ' J 2k 

1 1 «; 

x <D° : q i\<bt i] \+Z2m% 2 4-Z2/2)\p :q,v). 



X 



(17) 



From now on we will omit the color and spinor labels whenever this does not lead to ambiguities 
and replace the field-operator arguments z\ and z,2 by their non-vanishing components z\ and 
z 2 , respectively. Furthermore, if one uses k^ = x\p + and k 2 = X2Q~ to rewrite the and k 2 
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integrations in the amplitude (17) as integrations over the longitudinal momentum fractions x\ 
and x 2 , respectively, one arrives at, 



M MV =J dx,p + J ^e^ + ^ J dx 2 q- J i^e*^ 

x (D» :p'\V c (-zY/2)&M(+zY/2)\p : p,fi) H (x,p + ,x 2 q~) (18) 
x (D° : 9 '|^M t (+z+/2)^(-z+/2)|p : 4, v) . 

As in [1] for p — )■ A+ (/? — >■ A c ), the p — >■ D° (/? — y D°) transition matrix element is expected 
to exhibit a pronounced peak w.r.t. the momentum fraction. The position of the peak is approx- 
imately at 

in 

x =-^ = 0.68. (19) 
M 

From Eq. (9) one then infers that the relevant average momentum fractions jti and x 2 have to 
be larger than the skewness This means that the convolution integrals in (18) have to be 
performed only from ^ to 1 and not from to 1 . 

In the following section we will analyze the soft hadronic matrix elements in some more 
detail. 



4 Hadronic Transition Matrix Elements 

Compared to (10) the Fourier transforms of the hadronic matrix elements for the p — y D° and 
p — y D° transitions are rendered to Fourier integrals solely over and z 2 , respectively. Hence 
we have to study the integral 

p + J ^e m ? +z ^DV:p'\¥ c (-z^/2)<t> s W(+zi/2)\p : p,\x) (20) 

over the p — y D° transition matrix element and the integral 

q~ J ^e^ z 2(D° : 9 '|^ 5 M t (+z+/2)^(-z+/2)|p : q,v) (21) 

over the p — y D° transition matrix element instead of (11) and (12), respectively. 

We will first concentrate on the p — y D° transition (20) and investigate the product of field 
operators x ¥ c (—z~[/2)<i> s \- uc ^(+z~[/2). For this purpose we consider the c— quark field operator 
*P C in the hadron frame of the outgoing D°, cf. e.g. [10, 11], where the D° has no transverse 
momentum component. It can be reached from our symmetric CMS by a transverse boost 
[12, 13] with the boost parameters 

b + = (l-^)p+ and b ± = ^. (22) 
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In this hadron-out frame we write the field operator in terms of its "good" and "bad" LC com- 
ponents, 



1 



¥ c = ^(Y~Y + + Y + Y~) X V C = y ¥ c + + x ¥ c . 



(23) 



by means of the "good" and "bad" projection operators £P + = jYY + an d = jY + Y~, 
respectively. After doing that we eliminate the y~ appearing in ^ + and by using the 
c— quark energy projector 



£v(fc , 1 ,A 1 , )v(fc , 1 ,A 1 / ) =k[-Y-m c . 
K 

In the hadron-out frame it explicitely takes on the form 



(24) 



Zv{V l X)v(V l X)=k\ + Y- + ^ T Y + 
X', LK \ 



m 



C ) 



(25) 



since there the c— quark momentum is 



,J+ m c 



With those replacements the c— quark field operator becomes 

^ = ^Z{v(k[A{)(HKX)Y + v c ) 

1 A, 



+ 7 



v(F 1 ,A / )(v(^ 1 ,A 1 / )^) + 2m c ^]}. 



(26) 



(27) 



As in the case of pp — > A+A,, in Ref. [1] one can argue that the contribution coming from 
(v(k\,X[)Y +x ¥ c ) dominates over the one in the square brackets and thus the latter one can be 
neglected. Since this dominant contribution can be considered as a plus component of a four- 
vector, one can immediately boost back to our symmetric CMS where it then still holds that 



VJ/C(_Zl/2) = 2^ I v ^^ {WiM^-zi/l)) • 



(28) 



1 A,' 



Furthermore, one can even show that in (v(k[, X[)y +x ¥ c {-z^ /!)) on the r.h.s. of (28) only the 
"good" component of *P C (— z[ jT) is projected out, since 



v(k\, A 1 / )7 + ^'(-^/2) = v(k[, X[)y + ^ c {-z-J2) 



(29) 



Finally we note that such manipulations are not necessary for the scalar field operator 4> 5 M of 
the S[ud]— diquark. 
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p + J 


fd Zl 
2% 


e ix\p 


p + 
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X[ J 
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Putting everything together gives for the p — >■ D° transition matrix element (20) 
p + j ^-e^? + ^(D° : p'\*¥ c {-z[ /2)& S W (+^/2) \p : p,jU> = 

_ e mP + zi {]5 o . p / jv A /) (_ zr /2))*sM( +zr /2)|p : p,/i> . 

(30) 

Proceeding in an analogous way for the p — >■ D° transition matrix element, where the role 
of the + and — components are interchanged, we get for (21) 

q~ J ^e l ^ z 2(D° : /2)W{-4/2)\p : q, v) = 

|^ E / ^ 2 ^>° : ?'|^ t (+4/2) (r + (-z+/2)r. (k' 2 X)) u (k' 2 X) \P : v> . 

(31) 

Also here only the good components of the quark field are projected out on the r.h.s.. 

Using now (30) and (31) and attaching the spinors v(k[, X[) and U(k' 27 X 2 ) to the hard sub- 
process amplitude H by introducing 

H X[,l>, C*l 1*2) ■■= u (k'^X^ H (xip + ,x 2 q-) v (k' u X[) , (32) 
we get for the pp — > D°D° amplitude (18), 

M MV =^y 2 jdx, Jdx 2 H X[ ^ (* ,X 2 ) ^ ^ 

x v(*i, aot + p + / |^^ +zr : /in(-^r/2)* 5M (+^r/2)iP = p,m> 

(33) 

Introducing the abbreviations 

J*** := v(M, A0r + P + J %^ P+Z ~ l (D° ■■ P'\m-zi/2)^ s[ud] (+zi/2)\p : p,/i> (34) 
and 

: = / : q'\* S[ud] \+zl/2)%(-zl/2)\p : 4, v) y~u (k' 2 X) (35) 

for the pertinent projections of the hadronic transition matrix elements, we can write the hadronic 
scattering amplitude in a more compact form, 

AV = £ j dx x j dx 2 H X[ <K (x x ,x 2 ) ^— ^ jq^ . (36) 
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5 Overlap Representation of J^f^ 

In the following section we will derive a representation for the hadronic p — >■ D° and p — > D° 
transition matrix elements as an overlap of hadronic light cone-wave functions (LCWFs) for 
the valence Fock components of p and D° [6]. Since we only need them for x > i.e. in the 
DGLAP region, the hadronic transition matrix elements admit such a representation. For doing 
that we will make use of the Fock expansion of the hadron states and the Fourier decomposition 
of the partonic field operators in light cone QFT. 

At a given light cone time, say z + = 0, the "good" independent dynamical field components 
<§,s[ud\ an( j vj/<^ (— zj~ /2) of the S[ud\ — diquark and the c— quark, respectively, have the Fourier 
decomposition 



a (S[ud\ : k^, ki±) e 



-ik 



1 T 



+tf (S[ud] : k 1± ) e 



+ik 



1 2 



(37) 



and 



1 Ai - 



(38) 



The spinors u + and v + are the "good" components of the (anti-)quark spinors u and v, i.e., 
u + = £P + u and v + = t^? + v. The operators a and are the annihilator of an S [ud] -diquark 
and the creator of an 5 [ud] -diquark, respectively. The operator c annihilates a c-quark and the 
operator J t creates a c-quark. Their action on the vacuum gives the single-parton states 



a f (S[ud] : k+,k lA _) |0) 
(S[ud] : ^,k 1± ) |0) 
c f (c:^ + ,k' 1± ,AO|0) 
d f (c:*i+k' 1± ,^) |0> 



= \S[ud\:k+,k 1± ), 

= | %J] : &+ k 1± ) , 

= \c: £'+, k^), 

= |c:^+,k' 1± ,A0, 



(39) 
(40) 
(41) 
(42) 



which are normalized as follows 

(•••: £ /+ ,k' ± ,A'| •••:£ + ,k ± ,A) = 167T 3 fc+ 5 -Jfc+) 5 (2) (k' ± -k ± ) 5 A , >A . (43) 

In the case of S[ud] states no and no dx'.x appear. This normalization is in accordance with 
the (anti-)commutation relations 

a (S[ud] : k'+, k' ± ) , a f (S[ud] : k+, k ± )l = [ft : k' ± ) , £ f (S[i«/] : kj 

= 167r 3 fc+5(yt /+ -fc+)5^ (k' ± -k ± ) 



(44) 
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and 

{c (c : k ,+ , k' ± , A') , c t (c : *+ k ± , A)} = {d (c : k' ± , A') , ^ (c : *+ k ± , A)} 

= 16^ 3 fc+ 5 -*+) 8® (k' ± -k ± ) 5^ . 

In the Fock state decomposition hadrons on the light front are replaced by a superposition 
of parton states. Taking only into account the valence Fock state, the proton- and the D°-state 
in our quark-diquark picture are represented as 

\ P : P ,»)=Jd^ Wp (*,* L )^=L= (46) 
x I s [ud] ■ *P + , k_L +xp±) \u: (1 -x)p + , -k ± + (1 -x)p±, fi) 



and 



D° :/>=/ djc' Vd k'J - 1 -1= V (2A ; ) 
^ 7 167T 3 ^ V ±J Jx>(\-X>) V2ji K ' 



^x>(\-x>) v^tT (47) 
x I c : x'p' + , k' ± +x / p' ± , A') | w : (1 -iV + , ~ k l + (1 -*V±> - A '> > 
respectively, with normalization 
(• • • : p' ± (, A') | p + , p ± (, A)) = 16* V 6 {p' + - p + ) (p' ± -p ± ) (5 A , >A ) . (48) 

Also here, in the case of the pseudoscalar D° and D° states no A W and no 5 A / ;A appear. ty/ p and 

y/b are the LCWFs of the proton and the D°, respectively, which will be specified in Sec. 7. 
The LCWFs do not depend on the total momentum of the hadron, but only on the momentum 
coordinates of the partons relative to the hadron momentum. Those relative momenta are most 
easily identified in the hadron-frame of the parent hadron. Here we have assumed that the 
partons inside the proton and the D° have zero orbital angular momentum. The arguments of 
the light cone- wave functions are related to the average momenta and momentum fractions by: 

*i = TTT' kl± = kl± -TH—> (49) 

X2 r j x 2 A ± 

x 2 = yq^F , k 2± = k 2± + 1 + 2 , (50) 



fcl_L 


= k l± - 


k 2 ± 


= k 2 ± + 




= h± + 


k 2 ±_ 


= k 2A _- 



*i = Trf k' 1± = k l±+ 1 T ^^, (5i) 

v __ x 2 v - x 2 A ± 

x 2 — | _ £ ' K 2± _ K 2±. j ^ • P Z J 

Using the expressions above, we can write the hadronic matrix elements appearing in (33) as 

dxd 2 k 



^=-V2^J d ^^ WD {x\x^)^^m (53) 
x y^(x(x,£),k ± (k ± ,x,<f;)) 5(^i-^)5_ A;;M 
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for the p — > D° transition and 



dyd 2 l± — % 



VD(f(y,!;),V ± Q±,y,S)) 



I6n 3 Vy + t, 

x V P {y(y,Z),h(h,y,Z)) S(x 2 -y)S_^ v 

for the p — >■ D° transition. Here we have used that 

v (k[ , X[) r+v (k[ , -/i) = 2&'+ and u - v) y~w (4, ^2) = 2*£ - 
Collecting all pieces we finally get 

jUV f _ f _ __ 1 1 

M MV =— / rfxi / dx 2 H-n.- v (xi,x 2 ) 



(54) 



(55) 



X / 16^ Vto(^(^i»0» k' ± (k ± , xi, £))yp(*(*i,£)> k ± (k ± , xi, §)) 

/d 2 l 
^3 yfo (/(x 2 , £), ^2 ,$))V P {y(x 2 , £), 1±(I± ,*2 • 



(56) 



Furthermore, we can take advantage of the expected shape of the p — >■ D° (p — >■ Z)°) transi- 
tion matrix elements. Due to their pronounced peak around xo only kinematical regions in the 
hard scattering amplitude close to the peak position are enhanced by the hadronic transition ma- 
trix elements. For the hard partonic subprocess we therefore apply a "peaking approximation" 
i.e. we replace the momentum fractions appearing in the hard-scattering amplitude by xq. Then 
the hard scattering amplitude can be pulled out of the convolution integral and the pp — > D°D° 
amplitude simplifies further to 



Afyv =^2~#-m,-v(>o,*o) 



i dx v^hF / 



d 2 k 



167T 3 



(57) 



Yd (x (x, £ ) , k' ± (k ± ,x , £ )) x\f p (x{x, £ ) , k ± (k ± ,x , £ )) 
where the term in square bracket can be considered as a sort of generalized form factor. 



6 Hard Scattering Subprocess 

Before we start to specify the LCWFs occuring in this overlap representation of the p — >■ 
D° (p — > D°) transition, we will first calculate scattering amplitudes of the hard partonic 
S[ud] S[ud] — > cc subprocess within the peaking approximation. 

The hard- scattering amplitudes for the hard partonic subprocess, as shown in Fig. 3, is given 

by 

H X >.X> = 4 ("'ft* ( k 2^) r M v (k\M)) —^-2 ((-igsF s ) (*! -k 2 ) v ) . (58) 
y {ki+k 2 ) 
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immmmm 



S[ud] (fci) 



Figure 3: The hard scattering process on the partonic level S[wd]S[w<i] — > cc. 



4/9 is the color factor which we have attached to the hard- scattering amplitude. g s = y/4noc s is 
the "usual" strong coupling constant and F s denotes the diquark form factor at the gluon-diquark 
vertex. This diquark form factor takes care of the composite nature of the S[u d] -diquark and the 
fact that for large s the diquark should dissolve into quarks. We have taken the phenomenolog- 
ical form factor from Ref. [14], namely 



It is just the analytic continuation of a space-like form factor to the time-like region. The original 
space-like form factor was introduced in Ref. [15] (where it has been obtained from fits to the 
structure functions of deep inelastic lepton-hadron scattering, to the electromagnetic proton 
form factor and elastic proton-proton data at large momentum transfer). It should be remarked 
here that such a continuation is not unique, the form factor can acquire unknown phases when 
doing the continuation. But fortunately such phases are irrelevant w.r.t. the physics, which is 
the reason for taking the absolute value in (59). 

With the help of the peaking approximation we can express the subprocess amplitudes in 
terms of the kinematical variables of the full process. For the different helicity combinations 



Ql 



<2^ = 3.22GeV 2 , s>Qq. 



(59) 
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we explicitely have, 

. 2 , . 2 , 4 2M 
H++ = +4na s (x s) F s (x s) - -j= cos , 

= -4^5(^5) F 5 (jc§,s) ^ sin , 
= -4^a,(^5) F a (j^s) ^ sin , 
// = -47C(X s (xqs) F s (xqs) ^ ^cos . (60) 

7 Modelling the Hadronic Transition Matrix Elements 

In order to make numerical predictions we have to specify the LCWFs for the proton and the 
D°. We will use wave functions of the form 

_ a 2y k jW 

L ' * , (61) 

which can be traced back to a harmonic oscillator ansatz [17] that is transformed to the light 
cone [18]. In [19] it has been adapted to the case of baryons within a quark-diquark picture. 
According to [1,16] we write the wave functions of a proton in an S[ud]u Fock state as 

2 k2 

¥ P (x, k ± ) = N p xe~ a "^) (62) 
and the one of a pseudoscalar D° meson in a uc Fock state as 

..2 M 2 (x-*o) 2 -2 k l 

Vfo(*,kj_) =iV D e ^ g fl oi(T^y_ (63) 

Here x is the momentum fraction of the active constituent, the S[«<i]-diquark or the c— quark, 
respectively. The mass exponential in (63) generates the expected pronounced peak at x ^ xq 
and is a slightly modified version of the one given in [19]. 

In each of the wave functions, (62) and (63), we have two free parameters: on the one hand 
the transverse size parameter a p / D and, on the other hand, the normalization constant N p / D . The 
parameters can be associated with the mean intrinsic transverse momentum squared (k]_) p m of 
the active constituent inside its parent hadron and with the probability to find the hadron in 
the specific Fock state (or with the decay constant f p j D of the corresponding hadron). The 
probabilities and the intrinsic transverse momenta for the valence Fock states as given in (46) 
and (47) can be calculated as 

d 2 k\ 



Pp/D = J dxj I \]/ p/D (x, kj_) | 2 (64) 
and 

( k l) P /D = ^J dx J i6^ k i I Vp/d{x, kO | 2 , (65) 
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respectively. Inserting the wave functions (62) and (63) into (64) and (65) we obtain 



and 



Pd = 



Ni 



D 



32K 2 a 2 D 



N 2 
6A0n 2 a 2 p 



I n {a 2 D ) 



2la 2 p 



/t 2\ 1 hijal) 
2a D I n {a D ) 



(66) 



(67) 



where we have introduced the abbreviation 



.(4) := /' 
Jo 



dxx n (\-x) m exp 



-2alM^ X - X0) 



x(l — x) 



For the proton we use the same parameters as in [1] and [16]. We choose a p = 
for the oscillator parameter and P p = 0.5 for the valence Fock state probability. 



(68) 

l.lGeV 1 
Choosing 

P p = 0.5 for the proton may appear rather large at first sight. As a bound state of two quarks a 
diquark embodies also gluons and sea quarks and thus effectively incorporates also higher Fock 
states. Therefore a larger probability than one would expect for a 3 -quark valence Fock state 
and a larger transverse size of the quark-diquark state appear plausible. Choosing the parameter 
values as stated above we get for the proton 



280MeV and N n = 61.818 GeV 



-2 



(69) 



For the D-meson we fix the two parameters such that we get certain values for the valence 
Fock state probability Pq and the decay constant fy. The decay constant fo is defined by the 
relation 

(0 | ^"(0)/75^ c (0) | D° : p) = if Dl P . (70) 
Taking the plus component and inserting the fields as given in Sec. 5 we get (omitting phases) 



r f d 2 k ± 
2V6 / dx j^ Wd(x, fcjj = f D . 



such that 



\6K 2 a 2 D f D 



N D = 



(71) 



(72) 



2v / 6/n(fl| ) /2)' 

As value for the decay constant we take the experimental value fo = 206 MeV from Ref. [8], 
for the valence Fock state probability we choose P D = 0.9. This amounts to a D = 0.864GeV _1 . 
As values for the normalization constant and for the root mean square of the intrinsic transverse 
momentum of the active quark we then get 



\J (k^o = 383 MeV and No = 55.202 GeV~ 2 , 



(73) 



respectively. 
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Let us now turn to the issue of the error assesment w.r.t. the parameters. For the decay 
constant of the D°-meson we take = 206 ± 8.9 MeV as stated in Ref. [8]. The valence Fock 
state probability of the D-meson Po is varied between 0.8 and 1. We do not take into account 
the uncertainties of the parameters appearing in the proton LCWF They are small compared to 
the ones of the D-meson LCWF since they have been determined from detailed studies of other 
processes. The influence of the parameter uncertainties on the cross sections are indicated by 
grey error bands in Figs. 5 and 6. 

We now turn to the wave-function overlap as derived in Sec. 5. When taking the model wave 
functions (62) and (63) we explicitely get 

d 2 ki 



^N p N D (x + ^){x 2 -^)(l-x) 



1 



167T 2 

x exp 



(1 + 



2 M 2 (*-£-*0(l-£)) 2 
D (,-0(1-,) 



a 2 D (l-^ 2 (x + ^+a 2 (l+^ 2 (x-^ 

a 2 D a 2 p (l-x) 



exp 



a 2 D (l-^(x + ^ + a 2 M+^(x-^ 



(74) 



2.0 r 



2.0 r 




Figure 4: The wave-function overlap of Eq. (74) versus x at CMS scattering angle = (left 
figure) and 9 = x/2 (right figure). We show it for Mandelstam s = 30, 20 and 15 GeV 2 (solid, 
dashed and dotted curves). 

In Fig. 4 we show the wave-function overlap of Eq. (74) versus the momentum fraction x 
with the parameters chosen as stated above. First we observe that it is centered at x xq for 
vanishing CMS scattering angle. Next let us compare the left and the right panel. We see that the 
magnitude of the wave-function overlap is strongly decreasing with increasing CMS scattering 
angle 6. The wave-function overlap is also more pronounced in magnitude and shape in forward 
direction. Furthermore, when comparing the overlap for different values of Mandelstam s, we 
observe that in the more important forward scattering hemisphere the overlap is increasing in 
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magnitude with increasing CMS energy s, whereas at large scattering angles this behavior is 
reversed. 



8 Cross Sections 

The differential cross section for pp — > D°D° reads 



da P p^pOD° 1 da pp^nPnP 1 IAm , 
dQ. " 4^ SAMAm dt " 647? 7 A^" (75) 

where we have introduced 

4 



a o := iE I AVI 2 • ( 7 6) 



do-/dt [nb/GeV 2 ] do-/dt [nb/GeV 2 ] 




Figure 5: The differential cross section do p -^j^ o0 /dt versus | t' |. On the upper left, upper 
right and lower panel we show it for Mandelstam s = 15, 20 and 30 GeV 2 , respectively. 

In Fig. 5 the differential cross section dO p .^j^ oQ /dt is plotted versus 1 1' |, again for Man- 
delstam s = 15, 20 and 30 GeV 2 . The decrease of the cross section with increasing | t' \ can 
mainly be attributed to the wave-function overlap which gives rise to a generalized form factor 
(cf. Eq. (57)). This form factor enters the differential cross section to the fourth power. The 
forward direction is dominated by those amplitudes in which the helicities of the proton and an- 
tiproton (and also of the c and c quark) are equal. They go with cos 0. With increasing scattering 
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angle they compete with those in which proton and antiproton (and also c and c) have opposite 
helicities. The latter go with sin and dominate at 90°. If one looks at the energy dependence 

one observes that M ++ and M are suppressed by a factor 2M / y/s as compared to M + _ and 

M_ + (cf. Eq. (60)). This suppression is the reason that the cross section for pp — y D°D° is 
smaller than the one for pp — y A+A c in Ref. [1]. In the pp — y A+A c case the factor Mj yjs 
comes with those amplitudes which vanish in forward direction. When comparing the different 
panels of Fig. 5 one sees that the effect of the increase of the differential cross section with 
decreasing scattering angle becomes more pronounced for higher CMS energies. 

0.4 r 




Figure 6: The integrated cross section <^ p p_^jp versus Mandelstam s. 

In Fig. 6 we show the integrated cross section o versus Mandelstam s. It is of the or- 
der of 10 *nb, which is one order of magnitude lower than the integrated cross section for 
pp — y A+A c in Ref. [1]. This finding is in accordance with the diquark-model calculation of 
Ref. [16]. According to Ref. [16] larger cross sections (of the order of the pp —y A+A c cross 
section) are to be expected for the pp —y D + D reaction. This, however, requires to extend our 
handbag approach by including vector diquarks and will be the topic of future investigations. 
Our estimated cross section is about two orders of magnitude smaller than the predictions given 
in Refs. [20] and [21], where hadronic interaction models have been used. Whereas the authors 
of Ref. [20] determine their couplings of the initial proton to the intermediate and final charmed 
hadrons by means of QCD sum rules, the authors of Ref. [21] rather use SU(4) flavor symme- 
try. Though their predictions for the integrated pp —y D°D° cross section are comparable, they 
differ substantially in the pp — y D + D cross section which, in Ref. [20], is even smaller than 
our pp —y D°D° cross section. Such big discrepancies reveal the high necessity of experimen- 
tal data which allow to decipher between different dynamical models. Such experiments could 
also help to pin down the charm-quark content of the proton sea. A considerably higher cross 
section within our approach could only be explained if the charm-quark content of the proton 
sea was not negligible. 
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Appendices 

A Kinematics 

The four-momentum transfer A can be written as (cf. (1), (2)) 

)-m 2 (l-^) 
2p+{l-&) 



M 2 (l + ^-m 2 (l-^)+^A 2 /2 ■ 



(77) 



Note that A + = — A since p' — p = q — q'. 

In order to find expressions for the sine and the cosine of the CMS scattering angle we 
write the absolute value of the three-momentum and the momentum component into z-direction 
of the incoming proton as 

T Am ' P3= ~2 



|p|= ^A m , P3 = ^/a^-A 2 /. (78) 



and that of the outgoing D° as 

| P 'I=^Am, \p' 3 \= £y/Ab-Ay S , (79) 

respectively. 

Note that we have chosen the coordinate system in such a way that the z-component of the in- 
coming proton momentum is always positive. But that of the outgoing D° can become negative 
at large scattering angles due to the unequal-mass kinematics. This change of sign occurs when 
A 2 reaches its maximal value 

Ai max =*AL (80) 
which follows directly from (79). Then the CMS scattering angle can be written as 

a ( \ , f Pi 

9 — arccos — + arccos —^r 

VIpI/ VIp'I 

/ |Aj_|\ f\A ± . , 

= arcsin — — + arcsin — — - for forward scattering (81) 



2 IpI/ \ 2 \v'\ 

|A±|\ . /|AJ 

— — - - arcsin — — - 

2|p|y V2|p1 



= arcsin I ] — arcsin ( I +7t for backward scattering. 
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Using (81) the sine and cosine of the CMS scattering angle turn out to be 

1 



sinO = 



and 



cos = 




--A + sign (p' 3 ) \ Al 



A M~ 



SI 



g n (^ 3 ) 



(82) 



(83) 



respectively, where sign (p' 3 ) takes care of the kinematical situation of forward or backward 
scattering. 

Now we are able to express several kinematical variables in a compact form. Starting from 
the definition (2) of the average hadron momentum p and using (78), (79), (82) and (83) its plus 
component can be written as 

P + = \{P + +P + )= ^((P0+P3)+(P'0+P'3)) 



= [l + ^A 2 - Al/s + sign^) sJk\ - Al/s] 

= \^l[ 2 + \/ A2 m +A 2 M + 2A m A M cose] . 



(84) 



Note that in our symmetric CMS q = p + . For the skewness parameter t, we get 

P + -P' + 



A 2 -Aj/s- sign(p 3 ) ^/A 2 M - AjA 
2 + y/Al-Al/s + sign (p> 3 ) ^Ai-Al/s 



(85) 



A 2 — A 2 



A 2 + A 2 , + 2A m A M cos 2 + J A 2 + A 2 , + 2A m A M cos G 



m~ M 



Note that, as a consequence of the unequal-mass kinematics, £ cannot become zero, which is 
different from, e.g., Compton scattering where t, would be equal to zero in such a symmetric 
frame. For p' 3 > 0, however, t, is fairly small in our case and tends to zero for s — > °°. 

Now let us further investigate the Mandelstam variables and write them in a more compact form 
with the help of (78), (79), (82) and (83). Mandelstam t can be written as 

A 2 2^ 



t = 



l-£ 2 l-£ 2 

2 



(l + ^)M 2 -(l-^)m 2 ] 



A- ± s 



A 2 m + A 2 ,-2sign(/ 3 ) JAl - A 2 /sJA 2 M - A\/s 



(86) 



A 2 +A^-2A m A M COS0 
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It cannot become zero for forward scattering but acquires the value 



to := t(Al = 0,p' 3 >0) = - S -(A m -A M )\ 



and for backward scattering 



h := r(Ai = 0,p' 3 <0) = --(A m +A M y 



(87) 



(88) 



It is furthermore convenient to introduce a "reduced" Mandelstam variable t' that vanishes for 
forward scattering, 



t-to 



A m A M - sign(p' 3 ) J A 2 -A 2 JsJA 2 m - A 2 Js 



(89) 



= --A W A M (1 - cose) . 



Also the transverse component of A can easily be written as a function of the sine and the cosine 
of the scattering angle using (82) and (83), 



Ai=s 



A^A^sin 2 ^ 



or solving (89) for A , one finds 



A m+ A li + 2A m A M cos e 



, sA m A M + 1' 



(90) 



A± = ~ f s/4(A m + A M ) 2 + t r 
If we define Mandelstam u for forward scattering in an analogous way 

u := k (A 2 = 0,p 3 > 0) = -^-(A m + A M ) 2 
and for backward scattering 

Ul ■= u{A\ =<V 3 <0) = -^(A m -A M ) 2 , 
the sine and the cosine of half the CMS scattering angle can be written compactly as 



sin 



cos 



1 — cos to — t 



(!)- 



2 sA m A M ' 

1+COS0 Ul — u 

2 sA m A M ' 



(91) 

(92) 

(93) 

(94) 
(95) 



respectively. 
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B Light Cone Spinors 

For our purposes we use the light cone spinors [22, 23]. They read 



u(p,t) 



v(p,t) 



1 1 



-1 1 



fp++m/V2\ 
P±/V2 
p + —mj\fl 
\ P±/V2 J 

( - P *JV2 \ 
p + —m/y/l 
P1/V2 
\-p+-m/V2j 



for p 3 > and 



u(pA) = 



1 sign(p l ) 
21/4 a/ — 



v(p,t) = 



1 ^(p 1 ) 

2V4 a /~o~ 



( P1/V2 \ 
p +m/v / 2 

\p - m/ 1/2/ 



/ p -mj\fl\ 
-P±/V2 
—p~ —m/y/2 



u(p,i) 



v(p,i) 



1 1 



-1 1 

2^7^ 



/ -P1/V2 \ 
p + +m/ y/l 
P1/V2 
\-p + + m/\fl) 

fp+-m/V2\ 

P±/V2 
p + +m/ \/2 

\ P±/V2 J 



(96) 



(97) 



u(p,i) = 



1 -sign (p l ) 

2 1 / 4 



v(p,i) = 



1 -sign(p l ) 
2 i/4 



-P±/V2 
p +m/Vl 

V P±/v^ / 

(98) 
/ PI/V2 \ 
p —m/\/2 

\p-+m/V2j 
(99) 



for t? 3 < 0. They satisfy the charge-conjugation relation 

v(p,X) = iy 2 u*(p,X) (100) 

and are normalized as 

ux'{p)ux{p)=2m8 xix and v A , (p) v A (p) = -2m5 A / A . (101) 



C TDAs 

Following Ref. [5] the p ^ D° transition matrix element can be decomposed at leading twist 
into the following covariant structures, 

=+ [ dz l 



j&cs . =p + J -^-e'^ + ^ (D° : p' \ <F + (-z^ /2) <t> s ^ (+z[ /2) \p:p,v) 



= y 5 u(p,ii)Vi{xi, £,t) + 
+ u(p,n)Vi(x h £,t) + 



4 



M + m 



y 5 u(p, ju)V 2 (xi, £,t) 



(102) 



4 

M + m 



u(p,n)V 2 (x h £,t) 
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where we have introduced the p — )■ D° TDAs Vi, V 2 , Vi and V 2 . When evaluating the — > 
amplitude the hadronic transition matrix element (102) appears within the spinor product 



cf. (33). Expressed in terms of TDAs we thus have 



(103) 



^ =v(*i, A() 7 + 75 u(p, n)Vi(x u Z,t) + v(k[, X[) Y + T7^Y5u(p, n)V 2 (x h $,t) 



4 



M + m 



+ v(k[ , X[) y + u(p, n) V X (x u £ t t) + v(k[ , X[) 7+ 



M + m 



1-e 



v(p', A 1 / )r + 75"(p, /i)Vi(*i, §,0 + v(p', A() y+— ^y 5 w(p, /i)V 2 (*i, 5,0 



+ v( P ', x[) y+u(p, n) vi (xut,t) + v( P ', x[) r+- 



M + m 



M + m 

u(p,n)V 2 (x h £,t) 



(104) 



after making the replacement k[ = x\p' in the v— spinor. Evaluating the various spinor products 
which appear in (104) by using the light cone spinors of App. B gives 



^cs . 4 P + /*i-£Aj 



M + mV ! + 2 



(V 2 (xi,Z,t)+V 2 (xi,Z,t)) 



M + mV l + £ 2 



(105) 



(106) 



and 



+ -2p + y / x l -W^ + ^ Vi(x h Z,t)-Vi(x h Z,t) 
It, m 



+ ■ 



(v 2 (*i,£,f)+?2(*i,£,0) 



(107) 



1 + M + m 

^ =-2p + v / iT r Iv / l+l[^^i^, 0+Vi(jci,|,0 
2^ m 



+ ■ 



(v 2 (*i,§,0-V2(*i,$,0) 



(108) 



1+^M+m 

The TDAs V\ , V\ , V 2 an d V 2 can now t> e expressed as linear combinations of J£^+ , , 

and J^s + For our overlap representation of the hadronic transition matrix elements we have 

j^f + = j^s_ = o (cf Eq (55)) This means that v 2 = v 2 = and 



Vi = 



Ap + ^^l^l + l 



,st^ — tyt — 1_ , 1 , 



(109) 
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Vi = 



; , <yt\ + <7t_ i . 



(110) 



We further have = - 



(cf. Eq. (53)), so that we finally get 



Vi = 



ocS 



and Vi = 0. 



(HI) 
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